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Abstract 

To  evaluate  reaction  rates  for  making  hydrogen  from  methanol,  kinetic  studies  of  methanol  decomposition,  methanol  steam  reforming,  the 
water  gas  shift  reaction,  and  CO  selective  oxidation  have  been  performed.  These  reactions  were  studied  in  a  microreactor  testing  unit  using 
a  commercial  Cu-Zn0/Al203  catalyst  for  the  first  three  reactions  and  Pt-Fe/y-alumina  catalyst  for  the  last  reaction.  The  activity  tests  were 
performed  between  120  and  325  °C  at  atmospheric  pressure  with  a  range  of  feed  rates  and  compositions. 

For  methanol  decomposition,  a  simplified  reaction  network  of  five  elementary  reactions  was  proposed  and  parameters  for  all  five  rate 
expressions  were  obtained  using  non-linear  least  squares  optimization,  numerical  integration  of  a  one-dimensional  PFR  model,  and  extensive 
experimental  data.  Similar  numerical  analysis  was  carried  out  to  obtain  the  rate  expressions  for  methanol  steam  reaction,  the  water  gas  shift 
reaction,  and  CO  selective  oxidation. 

Combining  the  three  reactors  with  several  heat  exchange  options,  an  integrated  methanol  reformer  system  was  designed  and  simulated  using 
MATLAB.  Using  this  simulation,  the  product  distribution,  the  effects  of  reactor  volume  and  temperature,  and  the  options  of  water  and  air 
injection  rates  were  studied.  Also,  a  series  of  optimization  tests  were  conducted  to  give  maximum  hydrogen  yield  and/or  maximum  economic 
profit. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  interest  in  hydrogen  production  for  fuel  cell  appli¬ 
cations  is  steadily  increasing  due  to  environmental  and  na¬ 
tional  security  concerns.  When  hydrogen  is  used  as  an  energy 
source,  it  produces  no  CO,  SO2,  NOx,  VOC’s  and  also  pro¬ 
duces  no  carbon  dioxide.  Of  course,  a  hydrocarbon  is  still 
needed  to  make  hydrogen  for  fuel  cells,  however  a  hydro¬ 
gen  fuel  processor  can  be  advantageous  because  of  its  en¬ 
ergy  efficiency  and  can  lower  greenhouse  gases  compared 
to  the  direct  combustion  of  hydrocarbons.  Therefore  the  use 
of  fuel  cells  for  transportation  and  electric  power  could  re¬ 
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duce  both  toxic  air  pollutants  and  green  house  gas  emis¬ 
sions. 

For  mobile  application  of  fuel  cells,  carrying  a  compressed 
gas  cylinder  as  the  hydrogen  source  seems  to  be  the  most 
common  choice  today.  However,  there  are  several  economical 
ways  to  make  hydrogen  from  hydrocarbons  such  as  gasoline, 
alcohols  and  natural  gas,  etc.  Of  the  many  possible  hydro¬ 
carbon  sources,  it  is  not  easy  to  choose  the  most  promising 
one  because  of  factors  such  as  safety,  economy  and  infra¬ 
structure.  Among  the  various  kinds  of  hydrocarbon  choices, 
the  advantages  of  high  energy  density,  easy  availability  and 
safe  handling/storage  are  now  making  methanol  one  of  the 
most  promising  sources  of  hydrogen  for  fuel  cell  vehicles  and 
small  stationary  power  sources. 

The  integrated  system  of  methanol  steam  reforming  to 
produce  hydrogen  usually  includes  four  reactions  in  three 
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reactors:  methanol  decomposition  (1)  and/or  methanol  steam 
reforming  (2),  water  gas  shift  reaction  (3),  and  CO  selective 
oxidation  (4): 

CH3OH  **  CO  +  2H2,  AH298°  =  90.64 kJ  mol-1  (1) 

CH3OH  +  H20  **  C02  +  3H2, 

AH298°  =  49.47  kJ  mol-1  (2) 

co  +  h2o  co2  +  h2, 

AH298°  =  -41. 17  kJ  mol-1  (3) 

CO  +  (l/2)02  C02,  AH298°  =  -28.3  kJ  mol"1  (4) 

The  first  two  reactions:  methanol  decomposition  and 
methanol  steam  reforming  on  a  commercial  Cu/Zn0/Al203 
catalyst  were  investigated  in  our  previous  paper  [1].  For 
those  reactions,  a  comprehensive  kinetic  study  has  been  re¬ 
ported  by  Peppley  et  al.  [2].  In  their  paper  they  derived  three 
rate  expressions,  one  for  methanol  decomposition,  one  for 
the  WGS  reaction  and  one  for  methanol  steam  reforming. 
They  derived  Langmuir-Hinshelwood  rate  laws  from  de¬ 
tailed  surface  mechanisms,  resulting  in  12  elementary  re¬ 
actions  and  two  different  active  sites.  Their  resulting  rate 
expressions  contained  numerous  parameters  and  are  overly 
complicated  for  our  purposes,  adding  computation  time  with 
no  gain  in  predictability.  More  recently  there  have  been 
two  other  papers,  which  evaluate  the  reaction  rate  by  an 
overall  Arrhenius  type  simplified  model  [3-6].  However, 
there  have  been  few  papers  that  addresses  both  methanol 
decomposition,  methanol  steam  reforming,  and  the  known 
by-products. 

In  our  previous  paper  [1],  an  eight  step  reaction  pathway 
for  methanol  decomposition  was  presented,  and  it  was  shown 
that  eight  compounds  are  formed  when  methanol  decomposes 
in  the  absence  or  scarcity  of  water.  To  design  an  efficient 
methanol  reforming  system,  it  is  necessary  to  know  the  exact 
rates  of  methanol  decomposition  and  by-product  formation. 
In  this  paper  we  proposed  that  five  elementary  reactions  are 
sufficient  to  evaluate  all  components  in  methanol  decompo¬ 
sition.  When  water  is  added  to  methanol  feed,  the  activity  of 
catalyst  increases  and  by-product  formation  decreases  signif¬ 
icantly.  Also,  H2  yield  increases  and  CO  yield  decreases  with 
increasing  ratio  of  water  to  methanol.  Therefore,  to  evaluate 
the  overall  methanol  reaction  rate  it  is  plausible  to  add  the 
methanol  steam  reaction  as  a  sixth  reaction  when  water  is 
present. 

The  water  gas  shift  reaction  is  a  critically  important 
reaction  to  shift  carbon  monoxide  and  water  to  hydro¬ 
gen  and  carbon  dioxide.  Although  numerous  studies  of 
the  reaction  kinetics  and  mechanism  for  this  reaction  have 
been  reported  during  the  past  decades,  there  are  still  dis¬ 
agreements  as  to  what  is  the  active  site  and  what  is  the 
reaction  mechanism  [7-16].  Analyzing  various  proposed 
mechanisms,  there  is  a  distinct  conflict  between  the  “ad¬ 


sorptive  mechanism”  (5)  and  “regenerative  mechanism” 
(6  +  7). 

H2O  +  CO  ->  CO(ads)  +  H20(ads) 

^  [hit*]  —>  C02(ads)  T  H-2(ads) 


■*  co2  +  h2 

(5) 

H20  +  Red 

-*  H2  +  0* 

(6) 

CO  +  0,  -* 

C02  +  Red 

(7) 

Although  many  rate  expressions  have  been  derived  based 
on  these  mechanisms,  an  empirical  rate  expression,  rco  = 
/cFcoFh2o(1  —  /?)can  successfully  predict  the  rate  of  the 
water  gas  shift  reaction  (where,  the  reversibility  factor,/?  = 

^h2^co2/^co^h2o^p)  [17,18]. 

Regardless  of  the  reformer  design  and  the  size  of  the  water 
gas  shift  reactor,  small  amounts  of  CO  exist  after  hydrocarbon 
reforming  and  after  the  water  gas  shift  reactor.  These  small 
amounts  of  CO,  typically  less  than  1  mol%,  must  be  removed 
to  prevent  poisoning  of  the  fuel  cell  electrodes.  Among  the 
various  methods  to  remove  CO  selectively,  catalytic  oxida¬ 
tion  is  considered  as  one  of  the  most  plausible  and  economical 
options.  In  the  CO  selective  oxidation  reaction  system,  in  ad¬ 
dition  to  the  CO  oxidation  reaction  (4),  the  H2  oxidation  also 
occurs (8) 

H2  +  (1/2)02  —  H20,  AH°298  =  -242  kJ  mol"1  (8) 

The  feed  stream  of  a  PROX  reactor  in  a  methanol  fuel  re¬ 
former  system  is  composed  of  H2,  O2,  CO,  CO2  and  H2O. 
Therefore,  in  addition  to  the  reactions  of  CO  and  H2  oxida¬ 
tion,  it  is  also  necessary  to  consider  the  water  gas  shift  reaction 
[19].  To  accurately  predict  the  concentration  of  all  gas  com¬ 
ponents  in  a  PROX  reactor,  all  three  reactions  (CO  oxidation, 
H2  oxidation  and  (reverse)  water  gas  shift  reaction)  must  be 
considered  simultaneously.  This  three  reaction  system  and 
the  exact  rate  expressions  for  each  reaction  are  important 
components  in  the  optimization  and  control  of  commercial 
fuel  reformers. 


2.  Experimental 

For  the  three  reactions:  methanol  decomposition, 
methanol  steam  reaction,  and  the  water  gas  shift  reaction,  the 
Cu/Zn0/Al203  catalyst,  a  commercial  catalyst  manufactured 
by  Sud-Chemie  (Catalyst  no.:  EX-2248)  was  used.  The  cata¬ 
lyst  used  for  CO  selective  oxidation  was  the  Selectoxo  cata¬ 
lyst  manufactured  by  Engelhard  (Pt-Fe/Al203).  All  catalysts 
were  ground  and  sieved  to  a  particle  diameter  of  200-250  p^m 
to  eliminate  internal  diffusion  resistance.  All  reaction  tests 
were  performed  in  a  standard  catalyst  performance  evalua¬ 
tion  unit.  A  stainless  steel  (or  glass)  tubular  reactor,  1/2  in. 
in  diameter  and  12  in.  long  was  used  for  all  reaction  tests. 
To  ensure  isothermal  conditions  along  the  bed  length,  a  split 
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Table  1 


Experimental  conditions  of  each  individual  reaction 


Conditions 

Methanol  decomposition/steam  reforming 

Water  gas  shift  reaction 

CO  selective  oxidation 

Catalyst 

Cu/Zn0/Al203 

Cu/Zn0/Al203 

Pt— Fe/Al203 

Manufacturer 

Sud-Chemie  Inc. 

Sud-Chemie  Inc. 

Engelhard 

Loading  amount 

l.Og 

1-Og 

0.5  g  +  5  g  inert 

Reaction  temperature  (°C) 

120-325 

120-275 

80-280 

Reaction  pressure 

Atmospheric 

Atmospheric 

Atmospheric 

Liquid  feed  rate(mlh_1) 

0.5-8 

0.5-8 

— 

Gas  feed  rate  (seem) 

— 

70-80 

160-180 

GHSV  (1  h-1) 

280-4400 

6000 

1000-20000 

Reactor 

SS  316/glass  tube 

Glass  tube 

SS  316 

tubular  furnace  was  used  and  the  temperature  of  the  catalyst 
bed  was  measured  directly  by  a  1/16  in.  J-type  thermocouple. 

The  reaction  tests  were  performed  at  temperatures  be¬ 
tween  120  and  325  °C.  For  methanol  reforming  experiments, 
methanol  feed  rate  was  controlled  precisely  by  a  syringe 
pump,  74900  Series  (Cole  Palmer),  from  0.5  to  8mlh_1, 
giving  a  volume  of  vaporized  methanol  at  STP  flow  of 
277-4427  cm3  h_1 .  In  cases  when  water  with  the  methanol  is 
added,  the  injection  rates  were  controlled  also  by  the  syringe 
pump  from  0.5  to  8mlh_1.  The  catalyst  load  was  between 
0.25  and  1.0  g  and  the  GHSV  at  reaction  temperature  was 
controlled  between  1000  and  10,000  h-1. 

For  the  water  gas  shift  experiments,  the  feed  gas  stream 
was  a  1 :2  mixture  of  CO  and  hydrogen  to  simulate  the  condi¬ 
tions  exiting  a  methanol  reforming  unit.  For  the  PROX  exper¬ 
iments,  the  feed  gas  stream  contained  62-72%  H2,  0. 5-5.0% 
02,  2-17%  N2,  0. 5-3.0%  CO,  and  20-24%  C02  (dry  ba¬ 
sis).  The  catalyst  load  was  1.0  g  for  the  WGS  and  0.5  g  for 
PROX  and  the  GHSV  at  reaction  temperature  was  controlled 
between  1000  and  20,000  h_1 .  All  the  reaction  runs  were  per¬ 
formed  under  atmospheric  pressure.  The  experimental  con¬ 
ditions  of  this  reaction  system  are  summarized  in  Table  1 . 

The  effluent  of  the  reactor  was  maintained  at  120  °C  with 
heating  tapes  to  avoid  liquid  condensation  and  connected  di¬ 
rectly  to  a  CARLE  Series  S  gas  chromatograph,  which  uses 
a  Hydrogen  Transfer  System  (Pd  membrane)  for  hydrogen 
analysis.  This  is  a  specially  designed  GC  with  dual  TCD 
using  two  different  carrier  gases:  helium  and  nitrogen.  He¬ 
lium  is  the  preferred  carrier  gas  for  all  components  except 
hydrogen  while  nitrogen  is  the  proper  carrier  gas  for  hydro¬ 
gen.  Two  columns:  Alltech  Chemisorb  107  (80-100  mesh, 
6  ft  x  1/8  in.)  and  Supelco  Carboxen  1000  (60-80  mesh, 
15  ft  x  1/8  in.)  were  connected  in  series  to  analyze  the  con¬ 
densible  and  light  gas  components.  Eight  components:  wa¬ 
ter,  methanol,  dimethyl  ether,  methyl  formate,  H2,  CO,  CH4, 
CO2  were  measured  during  each  test  run.  Material  balances 
on  carbon  were  calculated  to  verify  measurement  accuracy. 
For  CO  selective  oxidation,  the  GC  chamber  was  maintained 
at  a  constant  temperature  of  50  °C  to  separate  N2  and  O2 
properly.  Five  components  H2,  O2,  N2,  CO  and  CO2  were 
measured  during  each  test  run.  Material  balances  on  carbon 
were  calculated  to  verify  measurement  accuracy,  and  for  all 
runs  reported  here  were  within  3%  of  closure. 


3.  Results 

3.1.  Methanol  decomposition  and  methanol  steam 
reaction 

Before  running  the  steam-methanol  reaction,  tests 
of  methanol  decomposition  were  conducted  using  pure 
methanol  in  the  absence  of  water.  When  methanol  decom¬ 
poses  over  the  Cu/Zn0/Al203  catalyst  it  produces  hydro¬ 
gen  and  CO  as  main  components  and  several  by-products 
such  as  methyl  formate,  carbon  dioxide,  dimethyl  ether  and 
methane.  Fig.  1  shows  a  product  distribution  of  all  the  com¬ 
ponents  leaving  the  reactor  as  a  function  of  temperature  at 
a  space  velocity  of  4400  h-1.  Methyl  formate,  which  was 
a  main  by-product  when  methanol  is  not  fully  decomposed, 
reached  its  maximum  concentration  of  about  1 0  mol%  at  tem¬ 
peratures  between  285  and  300  °C.  The  other  by-products, 
CO2,  dimethyl  ether  and  methane  were  formed  between  1 
and  3%.  For  the  kinetic  study  of  methanol  decomposition,  43 
runs  were  made  at  five  different  space  velocities.  All  of  these 
data  points  are  shown  in  Fig.  2. 

When  water  is  added  to  the  feed,  methanol  decomposes 
more  rapidly  at  lower  temperatures  and  shows  a  different  pat¬ 
tern  of  product  distribution.  Fig.  3  shows  the  conversion  of 


Reaction  Temperature  (°C) 

Fig.  1.  Methanol  decomposition  in  the  absence  of  water:  mole  fraction  of 
the  products  at  the  reactor  outlet  (methanol,  hydrogen,  carbon  monoxide, 
methyl  formate,  carbon  dioxide,  dimethyl  ether,  methane,  catalyst  loading: 
1.0  g,  GHSV:  4400  h_1). 
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Fig.  2.  Conversion  of  methanol  with  various  space  velocities  in  the  absence 
of  water  as  a  function  of  reaction  temperature  (1.0  g  of  Cu/ZnO/A^Os  cat¬ 
alyst,  methanol  feed:  Laboratory  Grade,  GHSV:  280-4400  h_1). 


methanol  versus  reaction  temperature  when  water  is  added 
to  the  methanol  feed.  Complete  methanol  conversion  occurs 
near  250  °C  when  water  addition  is  more  than  30  wt.%  of  the 
feed  (43  mol%).  Another  effect  of  water  addition  is  the  reduc¬ 
tion  of  all  the  observed  by-products:  dimethyl  ether  (DME), 
methyl  formate  and  methane.  Methyl  formate,  the  main  by¬ 
product  for  methanol-only  decomposition  was  reduced  sig¬ 
nificantly  as  water  was  added  and  no  methyl  formate  was 
detected  when  the  feed  contains  greater  than  43  mol%  water. 
Also  there  was  no  DME  formation  when  the  feed  had  more 
than  24  mol%  water  and  no  methane  was  detected  after  a  feed 
of  8.6  mol%  of  water. 

3.2.  Water  gas  shift  reaction 

Fig.  4  shows  the  water  gas  shift  conversion  of  carbon 
monoxide  with  various  H2O/CO  ratios  in  the  temperature 
range  of  120-250  °C.  As  a  consequence  of  equilibrium,  the 


Fig.  3.  Methanol  conversion  in  methanol  steam  reaction  as  a  function  of  re¬ 
action  temperature  (water  content  in  feed  =  1 .8, 8.6, 24, 43, 64  mol%,  catalyst 
loading  =  l.Og,  GHSV  =  1100h-1). 


Reaction 
Temp  (C) 


250 

220 

190 

155 

*  *  *Q*  •  • 

120 

• •  A-  *  * 

175 

Feed  H20/CO  ratio 

Fig.  4.  Water  gas  shift  reaction:  CO  conversion  vs.  feed  H2O/CO  ratio  (re¬ 
action  temperature:  175-250 °C,  catalyst  loading:  l.Og,  pressure:  1  atm, 
GHSV:  6100  lh-1). 

conversion  of  CO  consistently  increases  with  the  H2O/CO 
ratio  at  constant  temperature.  For  example  a  1 : 1  molar  feed 
ratio  and  220  °C,  the  conversion  reaches  70%  while  the  equi¬ 
librium  conversion  is  calculated  to  be  87%. 

3.3.  CO  oxidation 

Fig.  5  shows  a  plot  of  CO  oxidation  conversion  at  various 
O2/CO  feed  ratios  as  a  function  of  reaction  temperature  for 
the  Selectoxo  catalyst.  As  shown  in  this  figure,  the  conversion 
of  CO  increases  at  high  O2/CO  ratios  and  low  CO  concen¬ 
trations.  When  the  O2/CO  ratio  was  2.0,  the  CO  conversion 
was  more  than  95%  even  under  200  °C.  In  the  three  cases 
where  the  O2/CO  mole  ratio  is  near  1.0,  lower  CO  concen¬ 
trations  gave  higher  CO  conversions.  All  data  points  in  this 
figure  were  used  to  determine  the  rate  expressions  of  three 
reactions  in  the  PROX  system. 


— • — 02/C0=2.9/2.8  — b — 02/00=5.1/2.5  — ♦ — 02/00=1.6/2.8 
02/00=1.0/1.2  -  --A--  -  02/00=1.2/0.6  - -a  --  02/00=0.6/0.6 


Temperature  (deg  C) 

Fig.  5.  Conversion  for  CO  selective  oxidation  with  various  O2/CO  ratio  as 
a  function  of  reaction  temperature:  0.5  g  of  Selectoxo  (Pt-Fe/y-alumina) 
catalyst,  flow  rate:  167  seem,  H2:  64-75  mol%,  O2/CO  in  mol%,  1  atm,  no 
water  addition. 
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2C11  yO  It 
(M  ethanol) 

(CH3)20 

y  <3> 

i<  1  > 

{DME) 

\ 

+ 

Ft  20 

2 CH  yO{a) 

{M  ethoxy) 

l<  2  > 

CO  +  h2o 

4  <  8  > 

2 CH  20{a) 

( Fo  rm  aldehyde) 

co2  +  h2 

i<  4  > 

+  H2 


+  H 


2 


2  CO  +  2  H  2 


<r 


<5> 


CH  yOCHO(a)  _<6>_>  cO,+CH 

( MethylFormate ) 

l<  7  > 


CH  yOCHO 
( M  ethyl  Fo  rm  a  te ) 


Fig.  6.  Schematic  diagram  of  simplified  eight  steps  reaction  mechanism  of 
methanol  decomposition. 


3.4.  Kinetic  studies  and  analysis 


3.4.1.  Methanol  decomposition 

The  reaction  mechanism  of  methanol  decomposition  and 
methanol  steam  reforming  was  studied  in  a  previous  paper 
[1],  however  a  complete  kinetic  model  was  not  presented 
there.  Based  on  the  reaction  mechanism  shown  in  Fig.  6  an 
eight  step  reaction  model  of  methanol  decomposition  can  be 
written  for  the  five  elementary  reactions.  The  rate  expressions 
assumed  for  each  reaction  are  listed  in  Table  2. 

The  constants  k\  to  k$  of  Table  2  are  assumed  to  be  func¬ 
tions  of  temperature  in  the  expression  of  A;  exp(£;/R7).  The 
10  unknown  parameters  were  found  using  the  following  steps. 
( 1 )  A  2  and  E2  were  determined  from  the  DME  and  MeOH  data 
alone,  using  a  material  balance  on  DME.  (2)  A4,  A5,  E4  and 
£5  were  then  determined  using  MF  and  methane  balances. 
(3)  A3,  £3  were  determined  from  the  WGS  experiments  in¬ 
dependently.  (4)  A 1  and  E\  were  determined  using  a  material 
balance  on  MeOH.  For  step  (1)  and  (2)  the  expression  in  Eq. 
(9)  was  minimized. 


minimize 


N 


F  ( pej  -  pc0‘ 


i= 1 


Table  3 


Activation  energy  and  frequency  factors  for  methanol  decomposition 


Reaction 

Activation 

energy 
(kJ  mol-1) 

Frequency  factor 
(molgCat_1  s-1  atm*) 

CH3OH^ 

CO  +  2H2 

91.8 

5854 

2CH3OH  — ; 

►  ch3och3+h2o 

88.9 

125 

CO  +  H2O- 

co2 + h2 

43.7 

17.973 

2CH3OH  —5 

►  CH3OCHO  +  2H2 

82.2 

1650 

CH3OCHO 

-*  co2 + ch4 

45.6 

1.174 

*  For  second  order  reaction,  the  unit  is  (molgcat  1  s  1  atm  2). 


For  steps  (3)  and  (4)  the  expression  in  Eq.  (10)  was  mini¬ 
mized. 


minimize 


N 


f  =  J2  (*e.i  -  XcjY 


i= 1 


(10) 


where,  F  or  /  is  the  objective  function  of  optimization  for 
minimizing,  i  the  number  of  data  points  (i  =  43  for  DME, 
MF,  and  MD,  i  =  30  for  WGS);  pQ  the  experimental  partial 
pressure  of  a  component  at  the  reactor  outlet  (atm);  pc  the 
calculated  partial  pressure  of  a  component  at  the  reactor  outlet 
(atm);  Xe  the  experimental  conversion  of  a  component;  X4  the 
calculated  conversion  of  a  component. 

The  exit  partial  pressure  of  each  component  was  deter¬ 
mined  by  numerically  integrating  a  one-dimensional  isother¬ 
mal  PFR  model.  For  example,  the  MeOH  balance  was  written 
as  follows: 


dftCH3OH 

d  W 


-E i\  (  PR?pco 

■M  exp  (  -^r  )  [  Pch3oh 


RT 


^p,md 


A2 exp  I  - -  ]  P} 


RT 


ch3oh 


A4  exp  I  - -  )  P} 


RT 


ch3oh 


(11) 


where,  /1CH3OH  is  the  molar  flow  rate  of  methanol,  W  the 
catalyst  weight  (g);  A/  the  frequency  factor;  Ei  the  activation 
energy;  Pj  the  partial  pressure  of  the  component  i  (atm);  7fp?mp 
the  equilibrium  constants  of  methanol  decomposition. 

MATLAB  subroutine  functions  “ODE23”  and  “Lsqnon- 
lin”  were  used  for  numerical  integration  and  optimization. 
A  summary  of  the  rate  parameters  is  shown  in  Table  3.  The 
quality  of  the  fit  for  the  overall  model  is  demonstrated  in 
Fig.  7  by  comparing  the  observed  and  calculated  rate  of  ex- 


Table  2 


Reactions  and  rate  expressions  for  methanol  decomposition  systems 


Reaction 

Formula 

Rate  law 

Methanol  decomposition 

CH3OH^>CO  +  2H2 

—n  =  ki(PcH3OH  —  P^Pco) 

DME  formation 

2CH3OH  — >  ch3och3  +  h2o 

— ri  =  ^2(Fch3oh)2 

WGS  reaction 

co+h2o*>co2  +  h2 

— O  =  k3(Pco^H2o  —  K~l  Ph2Pco2) 

MF  formation 

2CH3  OH  -*  CH3  OCHO  +  2H2 

— r\  =  k4(PcH3OH)2 

Methane  formation 

ch3ocho^co2  +  ch4 

~r5  =^5^MF 
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Fig.  7.  Comparison  of  experimental  data  with  calculated  value  for  the  rate  of 
methanol  exiting  reactor  (temperature:  175-250  °C,  catalyst  loading:  1.0  g, 
pressure:  1  atm). 

iting  methanol  for  all  experiments.  Although  some  data  is 
slightly  skewed,  the  overall  quality  is  quite  good. 

A  detailed  confirmation  of  the  accuracy  of  the  methanol 
decomposition  kinetic  model  is  accomplished  by  comparing 
the  calculated  and  experimental  exiting  molar  rates  of  all  the 
reactants  and  products.  Figs.  8  and  9  show  the  fit  of  DME  and 
MF  concentrations  by  the  model,  and  indicate  a  high  degree 
accuracy. 

3.4.2.  Water  gas  shift  reaction 

Fig.  10  shows  the  comparison  of  experimental  exit  flow 
rate  of  CO  with  the  calculated  values  obtained  from  kinetics 
in  the  water  gas  shift  reaction.  The  WGS  kinetics  were  found 
independent  of  the  methanol  decomposition  reaction.  An  em¬ 
pirical  rate  expression  was  used:  rco  =  kPcoPa2oO-  ~  P), 
where  ft  is  the  reversibility  factor  of  the  reaction  defined  as 
P  =  PCOl  PUl  /  Fh2o  Pco  Kp  and  Xp  is  the  equilibrium  con¬ 
stant.  The  values  of  A  2  and  E2  are  shown  in  Table  3. 


Experimental  DME  molar  rate  exiting  reactor  (mmol/Hr) 

Fig.  8.  Comparison  of  experimental  rate  with  calculated  value  for  (temper¬ 
ature:  20-325  °C,  catalyst  loading:  1.0  g,  pressure:  1  atm). 


Experimental  MF  molar  rate  exiting  reactor  (mol/Hr) 


Fig.  9.  Comparison  of  experimental  data  with  calculated  value  for  methyl 
formate  (temperature:  120-325  °C,  catalyst  loading:  1.0  g,  pressure:  1  atm). 

Attempting  to  predict  the  experimental  data  for  the 
methanol  steam  reforming  tests  using  only  the  five  reac¬ 
tions  of  Table  2,  which  excluded  the  methanol  steam  reaction 
(CH3OH  +  H2O  CO2  +  3H2)  was  not  successful.  However, 
when  a  sixth  reaction  with  the  rate  —r^  =  kc(Pcn3OHPn2o) 
is  added,  the  model  fit  of  the  data  greatly  improved.  There¬ 
fore,  a  sixth  reaction  was  added  to  simulate  the  reaction  of 
methanol  steam  reforming.  The  rate  expression  for  this  sixth 
reaction  was  found  to  be: 

/— 70358\ 

r6  =  2.67  x  lCrexp  f — — — j  (Pcu3ohPu2o) 

(£6  =  J  mol-1 ,  P  =  aim)  (12) 

3.4.3.  PROX  reaction 

For  the  CO  PROX  system,  a  reaction  model  in  which  three 
reactions  (CO  oxidation,  H2  oxidation  and  the  water  gas  shift 
reaction)  occur  simultaneously  was  chosen  and  the  empirical 
rate  expressions  derived  from  the  numerical  analysis  are  as 


Fig.  10.  Comparison  of  experimental  data  with  calculated  value  for  the  water 
gas  shift  reaction  (temperature:  175-250  °C,  catalyst  loading:  1 .0  g,  pressure: 
1  atm). 
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CO  Conversion  (%)  Experimental 


Fig.  11.  Comparison  of  experimental  data  with  calculated  values  for  CO 
selective  oxidation.  A  0.5  g  of  Selectoxo  (Pt-Fe/y-alumina)  catalyst,  flow 
rate:  167  seem,  FC:  62-72 mol%,  O2/CO:  0.5-2.0,  reaction  temperature: 
100-300  °C,  pressure:  1  atm,  no  water  addition. 


follows: 

9  /— 33092  \  0c  _01 

-n  =  3.528  x  10  exp  (  RJ  j  Po2pco 

(£i=Jmol_1,  P  =  atm)  (13) 


— V2  =  2.053  x  10  exp 


— 18742  \ 
RT  ) 


p0.5 

^02 


(£2  =  Jmol  1 ,  P  —  atm)  (14) 


HEl  HE2 


Fig.  12.  Process  flow  diagram  for  methanol  reformer  optimization  HEl, 
HE2,  and  HE3  are  the  heat  exchangers  and  Ml,  M2,  and  M3  are  mixing 
points. 


using  MATLAB.  The  simulation  assumed  one-dimensional 
PFR  models  and  used  the  subroutine  functions  “ODE23”  for 
integration. 

Fig.  13  shows  the  composition  profile  in  the  simulated 
methanol  steam  reformer  operated  at  a  point  to  produce 
enough  H2  for  a  1  kW  fuel  cell.  In  this  simulation,  the  feed 
rate  was  lOmolh-1  (methanol  +  water)  and  the  water  con¬ 
tent  was  23.8  mol%,  the  reactor  temperature  was  250  °C  and 
the  reactor  contained  250  g  of  Cu/Zn0/Al203  catalyst.  When 
comparing  the  simulation  result  with  the  corresponding  ex¬ 
perimental  data,  all  of  the  exiting  flow  rates  are  well  matched 
with  the  experimental  values 


— r3  =  4.402  x  103  exp 


-34104 


RT  V 
(£3  =  J  mol 


(  ^co^h2o  - 


Pcp2  Pu: 

K? 


-1 


P  =  atm)  (15) 


The  quality  of  the  fit  is  shown  in  Fig.  11  by  comparing  the 
observed  and  calculated  CO  conversion  and  selectivity  for 
all  the  experiments.  As  shown  in  the  figure,  the  calculated 
values  show  good  agreement  with  the  experiments. 


4.2.  The  effects  of  reformer  reactor  volume  and 
temperature 

Fig.  14  is  a  contour  plot  of  catalyst  weight  in  the  reformer, 
reformer  temperature,  and  H2  yield  per  unit  methanol  feed, 
for  the  three  reactor  integrated  system.  For  this  simulation 
the  process  variables  of  temperatures  and  catalyst  loadings 
of  WGS  and  PROX  reactor,  amount  of  water  addition,  and  the 


4.  Simulation  and  optimization 

4.1.  Combined  reactions  and  integrated  system 

The  integrated  system  of  methanol  steam  reforming  in¬ 
cludes  three  reactors:  reformer,  WGS  reactor,  and  PROX  re¬ 
actor.  In  this  reactor  system,  water  can  be  added  to  any  or  all 
three  reactors  and  a  variable  amount  of  air  must  be  added  to 
the  PROX  reactor.  Also  an  energy  balance  must  be  satisfied 
for  each  reactor  that  includes  heats  of  reaction,  latent  heats, 
sensible  heats,  and  heat  transfer.  Fig.  12  shows  the  simpli¬ 
fied  process  flow  diagram  of  the  methanol  reformer  system 
used  for  process  simulation  and  optimization.  In  this  inte¬ 
grated  system,  three  reactors,  three  mixing  points,  three  heat 
exchangers,  10  reactions,  and  10  components  of  reactants 
and  products  are  considered.  Based  on  the  reaction  kinet¬ 
ics  obtained  above,  a  process  simulation  program  was  coded 


Catalyst  Weight  (gr) 

Fig.  13.  Simulation  of  methanol  steam  reaction:  product  distribution 
through  the  reactor  bed  (feed  rate:  10  gmolh-1  of  methanol,  water  content: 
23.8  mol%,  inlet  temperature:  250  °C,  catalyst:  250  g,  isothermal). 
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Fig.  14.  Simulation  of  integrated  system  of  methanol  fuel  reformer:  effect  of 
reformer  volume  and  temperature  on  hydrogen  yield  (reactor  temperatures 
-  WGS:  220  °C,  PROX:  200  °C;  catalyst  loading  -  WGS:  500  g,  PROX: 
40  g;  water  addition  -MD:  10  molh-1,  WGS:  5  molh-1,  PROX:  3  molh-1, 
O2 /CO  ratio:  1.1). 

PROX  O2/CO  ratio  were  fixed  at  conditions  shown  in  Table  4. 
As  shown  in  Fig.  14,  the  hydrogen  yield  increases  expectedly 
with  increasing  reformer  catalyst  weight  and  increasing  re¬ 
former  temperature.  However,  the  CO  concentration  exiting 
the  integrated  system  (exit  of  PROX  reactor)  also  increases 
with  increasing  temperature  and  reformer  catalyst  volume, 
as  shown  in  Fig.  15.  To  meet  a  specific  CO  tolerance  speci¬ 
fication,  for  example  30  ppm,  it  is  expected  that  an  optimum 
condition  of  reformer  temperature  and  volume  may  exist. 

4.3.  The  effects  of  water  addition 

To  produce  hydrogen  from  methanol  or  hydrocarbon 
sources,  the  addition  of  water  is  always  advantageous.  But 
the  location  and  amount  of  water  to  add  is  not  as  obvious.  To 
add  water  to  the  reformer  system,  three  locations  are  possi¬ 
ble:  the  inlet  of  the  reformer  with  the  methanol  feed,  the  inlet 
of  the  WGS  reactor,  and  the  inlet  of  PROX  reactor.  Fig.  16 
is  a  contour  plot  of  hydrogen  yield  and  Fig.  17  is  the  exiting 
CO  concentration  versus  water  addition  rates  to  the  reformer 
and  the  WGS  reactor.  For  these  plots,  the  variables  of  re¬ 
actor  volume,  temperature,  and  O2/CO  ratio  are  fixed  at  the 

Table  4 


Base  conditions  of  the  variables  for  simulation  and  optimization 


Types 

Description 

Unit 

Conditions 

Reactor  size 

Methanol  reformer 

g  catalyst 

250 

Water  gas  shift  reactor 

g  catalyst 

500 

PROX  reactor 

g  catalyst 

40 

Reactor 

Methanol  reformer 

°C 

300 

Temperature 

Water  gas  shift  reactor 

°C 

220 

PROX  reactor 

°c 

200 

Water 

Methanol  reformer 

mol  h- 1 

10 

Injection 

Water  gas  shift  reactor 

mol  h- 1 

5 

PROX  reactor 

molh-1 

3 

02/CO  ratio 

PROX  reactor 

— 

1.1 

Temp  of  the  reformer  (deg  C) 


Fig.  15.  Simulation  of  integrated  system  of  methanol  fuel  reformer:  effect 
of  reformer  volume  and  temperature  on  CO  outlet  concentration  (reactor 
temperatures  -  WGS:  220  °C,  PROX:  200  °C;  catalyst  loadings  -  WGS: 
500 g,  PROX:  40 g;  water  addition  -  MD:  10 molh-1,  WGS:  5 molh-1, 
PROX:  3 molh-1,  02/CO  ratio:  1.1). 

conditions  in  Table  4.  These  figures  show  that  water  is  advan¬ 
tageous  for  both  hydrogen  yield  and  CO  outlet  concentration. 
Therefore,  the  proper  amount  of  water  should  be  determined 
after  considering  the  energy  penalty  imposed  by  the  water. 

4.4.  Reformer  optimization  and  economic  analysis 

Critical  to  the  system  design  is  finding  the  optimum  re¬ 
actor  sizes  and  operating  conditions  of  temperature,  amount 
of  water  added  to  each  reactor,  and  the  air  to  CO  ratio  to  the 
PROX  reactor.  For  example,  a  bigger  WGS  reactor  results  in 
a  smaller  PROX  reactor,  which  results  in  an  important  design 
trade  off  between  the  size  of  the  WGS  reactor  and  the  size  of 
the  selective  oxidation  reactor. 

To  make  a  compact  and  efficient  fuel  reformer,  many  vari¬ 
ables  must  be  optimized  simultaneously.  Based  on  the  kinet- 


Fig.  16.  Simulation  of  integrated  system  of  methanol  fuel  reformer:  effect 
of  water  addition  on  hydrogen  yield  (reactor  temperatures  -  MD:  300  °C, 
WGS:  220  °C,  PROX:  200  °C;  catalyst  loadings  -  MD:  250  g,  WGS:  500  g, 
PROX:  40  g;  water  addition  -  PROX:  3  molh-1,  02/CO  ratio:  1.1). 
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Fig.  17.  Simulation  of  integrated  system  of  methanol  fuel  reformer:  effect 
of  water  addition  on  product  CO  concentration  (reactor  temperatures  -  MD: 
300  °C,  WGS:  220°C,PROX:  200  °C;  catalyst  loadings -MD:  250  g,  WGS: 
500  g,  PROX:  40  g;  water  addition  -  PROX:  3  mol  h- 1 ,  O2/CO  ratio:  1.1). 

ics  of  the  reactions  and  our  simulation  program,  a  number 
of  optimization  tests  for  the  integrated  reformer  system  were 
performed.  For  this  optimization,  an  approximate  economic 
profit  was  used  as  the  objective  function.  Profit  is  defined  as 
the  difference  between  hydrogen  revenue  and  the  fixed  plus 
operating  cost  of  production.  In  this  optimization  the  tem¬ 
peratures  of  the  reactions  were  limited  from  100  to  350  °C 
to  avoid  condensing  water  and  thermal  degrading  of  the  cat¬ 
alysts.  Also  the  total  rates  of  water  addition  were  limited  to 
20  molh-1,  which  is  imposed  by  the  humidity  limits  of  the 
fuel  cell. 

For  hydrogen  revenue,  the  current  prices  of  cylinder  sup¬ 
plied  hydrogen  were  used.  The  prices  are  summarized  in 
Table  5.  Depending  on  purity,  the  hydrogen  price  varies  from 
0.05  to  1.12$  gmol-1 .  A  hydrogen  purity  factor  was  derived 
from  this  data  using  the  relation  between  H2  purity  and  prices. 
According  to  the  CO  outlet  concentration,  the  purity  factor 
(pufc)  is  expressed  as  follows: 

pufc  =  1.1983  exp(-0.0052CO  ppm)  (16) 

The  revenue  from  producing  hydrogen  was  calculated  by 
multiplying  the  hydrogen  price,  times  the  purity  factor,  times 
the  hours  of  reformer  operation  times  the  hydrogen  flow  rate. 

The  reactor  costs  were  determined  from  the  mass  of  cata¬ 
lyst  times  a  cost  factor  of  1  for  the  WGS  and  MSR  reactors 

Table  5 


Hydrogen  price 


Supplier 

Grade 

H2  purity 

Price  ($  kg  1 ) 

Scott  specialty  gases ' 

Research 

99.9999 

560 

Ultra-high  purity 

99.999 

150 

Pre-purified 

99.99 

84 

Airgas 

High  purity 

99.9 

44 

Industrial 

99 

25 

*  Based  on  cylinder  size  K:  262  ft3- 


and  5  for  the  PROX  reactor.  Energy  costs  for  heating  were 
calculated  assuming  methanol  is  combusted  to  provide  the 
necessary  heat  and  that  cooling  was  free.  These  are  opti¬ 
mistic,  but  a  place  to  start  our  analysis.  Factors  related  to 
hydrogen  price,  reactor  costs,  and  energy  costs  are  summa¬ 
rized  in  the  Table  6.  To  maximize  the  profit  function,  the 
optimization  program  minimized  (cost  —  revenue)  using  the 
“fmincon”  function  in  MATLAB  and  the  subroutines  previ¬ 
ously  discussed. 

The  methanol  decomposition  reaction  is  endothermic  and 
is  favored  at  high  temperatures  while  the  water  gas  shift  reac¬ 
tion  is  favored  at  low  temperatures.  However  the  reforming 
reactor  must  not  run  too  hot  to  prevent  catalyst  deactivation 
and  high  energy  costs.  Also  the  water  gas  shift  reactor  has 
an  optimum  temperature  because  the  rate  of  the  reaction  is 
kinetically  slow  at  low  temperatures  but  favored  thermody¬ 
namically.  CO  selective  oxidation  is  strongly  exothermic  with 
little  or  no  thermodynamic  limitations,  however  if  tempera¬ 
tures  are  too  high  the  reverse  water  gas  shift  reaction  rate 
increases  in  the  PROX  reactor  and  produces  CO  rather  than 
consuming  it. 

To  test  the  optimization  of  the  profit  function  it  was  easi¬ 
est  to  fix  several  of  the  10  operating  variables  in  Table  6  and 
search  for  optimum  values  for  the  remaining  variables.  The 
feed  rates  and  reactor  sizes  (except  the  reforming  reactor) 
were  fixed  at  those  values  in  Table  6  and  the  reactor  tem¬ 
peratures  that  gave  the  highest  profit  were  calculated  and  are 
plotted  versus  reformer  size  in  Fig.  18.  As  shown  in  this  fig¬ 
ure,  the  operating  temperature  of  the  reformer  is  optimum  in 
the  range  of  300-350  °C  while  the  optimum  temperatures  of 
WGS  and  PROX  reactor  are  very  constant  at  220  and  250  °C, 
respectively. 

Water  and  feed  rates  were  also  tested  for  optimum  be¬ 
havior.  In  Fig.  19,  optimum  water  addition  rates  were  ob¬ 
tained  to  maximize  profit  as  a  function  of  the  size  of  the 
methanol  steam  reformer.  All  other  values,  size  of  reactors, 
temperatures  of  reactors  and  O2/CO  in  PROX  were  fixed  at 
those  values  in  Table  6.  As  shown  in  Fig.  19,  there  exist  op- 


Reformer  -  -  -  WGS - PROX 


Catalyst  weight,  reformer  (gr) 


Fig.  18.  Optimum  operating  temperatures  of  the  reformer,  WGS  reactor, 
and  PROX  reactors  with  regard  to  the  size  of  the  reformer  (catalyst  load¬ 
ing  -  WGS:  300 g,  PROX:  50 g;  water  addition  -  MD:  10 mol h"1,  WGS: 
5  molh-1,  PROX:  3  molh-1,  O2/CO  ratio:  1.1). 
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Table  6 


Fixed  conditions  of  the  variables  in  the  optimization 


Symbol 

Description 

Unit 

Values 

h2price 

Hydrogen  price  (99.995) 

$  gmol-1 

0.25 

pufc 

Purity  factor 

— 

1.1983  exp(— 0.0052COppm) 

rxlcf 

Cost  for  methanol  steam  reformer 

$  Scat 

1 

rx2cf 

Cost  for  water  gas  shift  reactor 

$  Scat 

1 

rx3cf 

Cost  for  PROX  reactor 

$  Scat 

5 

qcl 

To  heat  feed  to  reformer  temperature 

$  °C_1  gmol-1 

1.41  x  1(T6* 

qc2 

To  supply  reaction  heat  for  MD 

$kJ-‘ 

1.72  x  10“5** 

qc3 

To  cool  reformate  to  WGS  temp. 

Free 

— 

qc4 

To  remove  reaction  heat  of  WGS 

Free 

— 

qc5 

To  cool  reformate  to  PROX  temp. 

Free 

— 

qc6 

To  remove  reaction  heat  of  PROX 

Free 

— 

*  Based  on  methanol  price,  0.64  $  gal  1  and  heat  of  combustion,  397  kJ  mol  1 
**  Based  on  the  heat  capacity  of  methanol  at  20  °C,  82.0  J  mol-1  K_1 . 


Reformer - WGS - PROX 


Catalyst  weight,  reformer  (gr) 


Fig.  19.  Optimized  water  addition  rates  to  reformer  WGS  reactor,  and  PROX 
reactors  with  regard  to  the  size  of  the  reformer  (catalyst  loading  -  WGS: 
300 g,  PROX:  50 g;  water  addition  -  MD:  10 mol h-1,  WGS:  5molh_1, 
PROX:  3  molh-1,  02 /CO  ratio:  1.1). 

timum  levels  of  water  addition  for  methanol  steam  reformer 
(13-15.5  gmolh-1)  and  WGS  reactor  (5.2-1 1.5  gmolh-1) 
while  a  lesser  effect  of  water  addition  exists  for  the  PROX 
reactor. 

5.  Conclusion 

The  reaction  kinetics  for  the  six  reactions  in  a  methanol 
steam  reformer,  the  one  reaction  in  the  water  gas  shift  re¬ 
actor,  and  the  three  reactions  in  a  CO  selective  oxidation 
reactor  were  measured  and  all  10  rate  expressions  were  ob¬ 
tained  from  a  large  set  of  experimental  data  from  commercial 
Cu/Zn0/Al203  and  Pt-Fe  catalysts.  Applying  these  empir¬ 
ical  relationships  in  simple  plug  flow  reactor  models,  it  is 
possible  to  gain  insights  into  the  design  and  operation  of  the 
reformer,  the  WGS,  and  the  PROX  reactors,  as  well  as  the 
integrated  system.  Using  these  rate  expressions  the  system 
was  simulated  and  optimized  to  maximize  profit. 

The  simulation  and  optimization  results  show  that  the  per¬ 
formance  of  the  integrated  system  is  greatly  affected  by  the 
size  of  the  reformer  and  not  sensitive  to  the  temperature  of  the 


WGS  reactor  or  PROX  reactor.  For  best  performance,  the  wa¬ 
ter  gas  shift  reactor  should  be  operated  in  the  range  of  220  °C 
regardless  of  other  process  conditions.  For  the  CO  oxida¬ 
tion  reactor,  the  operating  temperature  and  reactor  size  have 
less  impact  on  the  performance  of  the  reactor,  but  O2/CO  ra¬ 
tio  should  be  maintained  higher  than  stoichiometry  to  avoid 
high  CO  concentrations  in  the  final  product.  Although  our 
optimization  efforts  are  simplistic  and  brief,  their  usefulness 
is  significant.  The  results  of  this  study  are  expected  to  be  a 
critical  part  of  the  overall  design,  optimization  and  control  of 
fuel  reformer  systems. 
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